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Abstract. We present results of a population synthesis study aimed at examining the role of spin-kick alignment in producing 
a correlation between the spin period of the first-born neutron star and the orbital eccentricity of observed double neutron 
star binaries in the Galactic disk. We find spin-kick alignment to be compatible with the observed correlation, but not to 
alleviate the requirements for low kick velocities suggested in previous population synthesis studies. Our results furthermore 
suggest low- and high-eccentricity systems may form through two distinct formation channels distinguished by the presence 
or absence of a stable mass transfer phase before the formation of the second neutron star. The presence of highly eccentric 
systems in the observed sample of double neutron stars may furthermore support the notion that neutron stars accrete matter 
when moving through the envelope of a giant companion. 
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INTRODUCTION 

Recent observations of single and binary pulsars have 
sparked new questions and challenged accepted ideas 
on the formation of neutron stars (NSs) and the nature 
of supernova (SN) kicks. In particular, measurements of 
pulsar radio emission polarization and pulsar wind nebu- 
lae symmetry axis directions have provided increasingly 
compelling evidence for the alignment ofpulsar proper 
motions and pulsar rotation axes [e.g-QlEHIIlElla]- As- 
suming the proper motion and pulsar spin axis directions 
are representative of the natal kick and NS progenitor ro- 
tation axis, the alignment suggests that natal kicks are 
preferentially aligned with the progenitor's rotation axis. 
Moreover, the increasing sample of double neutron star 
(DNS) binaries has revealed a possible correlation be- 
tween the spin period P sp ; n of the first born NS and the bi- 
nary orbital eccentricity e (see Fig.[TJ |7, 8]. Such a rela- 
tion arises naturally for symmetric SN explosions, but is 
highly constraining for asymmetric explosions. Our aim 
in this paper is to examine the role of spin-kick alignment 
in establishing the observed P S p m -e correlation. 

DOUBLE NEUTRON STAR FORMATION 

According to our current understanding, DNSs form 
from binaries in which, possibly after some mass ex- 
change, both stars are massive enough to undergo core- 
collapse into a NS. After the formation of the first NS, the 
binary evolves through a high-mass X-ray binary phase 
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FIGURE 1. Orbital eccentricities and spin periods of ob- 
served DNSs in the Galactic disk. The dashed line represents 
the best-fitting log-linear curve e = — 1.7 + 1.21og,P S pj n . 



during which the NS accretes matter from the wind of 
its companion. When the companion exhausts its nuclear 
fuel, it swells up to giant dimensions and initiates a dy- 
namically unstable mass transfer phase. As the NS moves 
through the expanding envelope, orbital energy and an- 
gular momentum are transferred to the envelope, spin- 
ning it up until it is expelled from the system (assum- 
ing enough orbital energy is available for the expulsion). 
The binary emerges from this "common envelope" phase 
consisting of a NS and the helium core of its former gi- 
ant companion. If the NS is able to accrete matter while 
moving through the common envelope, the associated ac- 
cretion of angular momentum may spin it up to subsec- 



ond rotation periods. After the common-envelope phase, 
the giant's exposed helium core evolves further until it, 
in turn, explodes in a SN and forms a second NS. De- 
pending on the mass of the helium core and the size of 
the orbit right before the SN explosion, the formation of 
the second NS may be preceded by a dynamically sta- 
ble mass transfer phase during which accretion spins-up 
the first-born NS to millisecond periods. This formation 
channel is commonly referred to as the standard DNS 
formation channel 131 ■ 

The orbital characteristics of observed DNS binaries 
are determined by the mass loss and SN kick accom- 
panying the birth of the second NS, and the subsequent 
loss of orbital energy and angular momentum through 
gravitational radiation. In the case of an asymmetric 
SN explosion, the orbital semi-major axis A and eccen- 
tricity e right after the explosion depend on the mass 
loss and kick imparted to the NS through the conserva- 
tion laws of orbital energy and angular momentum [e.g. 
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Here G is the Newtonian gravitational constant, Mi and 
M 2 are the masses of the first- and second-born NS, Vo = 
[G(M\ +M ) /Aq] x I 2 is the relative orbital velocity of the 
second NS's progenitor right before its SN explosion, Mq 
is the pre-SN mass of the second NS's progenitor, Aq is 
the pre-SN orbital separation, and Vk is the magnitude 
of the natal kick velocity. The angles and define the 
direction of the natal kick velocity: 9 is the polar angle 
between Vk and Vo, and the corresponding azimuthal 
angle in the plane perpendicular to Vo (with = 7r/2 
corresponding to the direction from the first NS to the 
progenitor of the second NS). We note that in writing 
down Eqs. ©-©, the pre-SN orbit was assumed to be 
circular, as expected from the strong tidal forces acting 
during the common envelope occurring in the course of 
DNS formation. 

In the case of a symmetric SN explosion, no natal kick 
is imparted to the NS and the post-SN orbital elements 
are uniquely determined by the mass loss from the sys- 
tem. Equations 0}-© reduce to [e.g. [l5i[T6ll 
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FIGURE 2. Mass accreted by the NS as a function of the 
helium star mass at the end of a stable mass transfer phase in 
NS-helium star binaries (data from fl7l[T8l[T9ll ). The amount 
of mass accreted increases with decreasing donor star mass 
because the duration of the mass transfer phase increases with 
decreasing donor mass. The dashed line represents the best- 
fitting log-linear curve logAMj = —0.36 — 1.35Mo- 



THE P spin - e RELATION 

Faulkner et al. JUl noted that, in the absence of SN kicks, 
a relation between P sp i n and e arises naturally if the de- 
gree of spin-up of the first-born NS correlates with the 
mass of the progenitor of the second-born NS. Such a 
correlation is expected for binaries in which a stable mass 
transfer phase precedes the SN explosion of the progeni- 
tor of the second NS. The reason for this is twofold. First, 
less massive stars evolve more slowly than more mas- 
sive stars, yielding longer-lived mass-transfer phases and 
thus more time to accrete matter and spin the first NS 
up to shorter spin periods (see Fig.©. Second, in the ab- 
sence of SN kicks, the binary orbital eccentricity after the 
SN explosion forming the second NS is proportional to 
the amount of mass lost during the SN (see Eq.|4|. Less 
massive progenitor stars therefore not only imply more 
spin-up but also lower orbital eccentricities. Dewi et al. 
Il20ll have shown this theoretically expected correlation 
to persist if small kicks of the order of a few 10 km s _1 
are imparted to the second NS at birth. Larger kicks of a 
few lOOkms -1 , typical for isolated radio pulsars, spread 
out the range of post-SN orbital eccentricities, destroying 
any correlation between the first NS's spin period and the 
binary orbital eccentricity. 

The requirement of small kicks to reproduce the ob- 
served P sp i n -e correlation is puzzling considering the 
stringent lower limits on the kick velocities required to 
explain the observed binary properties and space veloc- 
ities of the relativistic DNS binaries PSRB 1534+12 and 
PSRB1913+16(see JHHH] and the paper by I. Stairs in 
this volume). Yet, these binaries follow the P %v \ a -e cor- 



relation shown in Fig.Q] A revision of the low kick sce- 
nario proposed to explain the P S pin-<? relation is therefore 
in order. 



POLAR KICKS 

In the absence of any theoretical models arguing other- 
wise, the formation of single and binary NSs has so far 
primarily been studied assuming all kick directions to 
be equally likely (some tentative considerations of non- 
isotropic kicks are presented in 112 ll 12311 1. However, in 
recent years, observational evidence that NS kicks may 
be preferentially aligned with NS rotation axes has be- 
come increasingly convincing: evidence for alignment 
between pulsar proper motion directions and pulsar wind 
nebula symmetry axes was presented by Wang et al. 
and Ng & Romani JH01, while Johnston et al. Qll^] and 
Rankin [5] presented evidence for alignment between 
pulsar proper motion directions and polarization vectors. 
Assuming the proper motion directions are representa- 
tive of the kick directions, these alignments suggest a 
correlation between pulsar rotation axes and kick direc- 
tions. In what follows, we furthermore assume the rota- 
tion axes of the pulsar and its progenitor to be parallel, 
so that the above evidence can also be considered as evi- 
dence for alignment between NS kick directions and NS 
progenitor rotation axes. We will refer to such kicks as 
"polar kicks". 

If SN kicks are restricted to be along the rotation axis 
of the NS progenitor {9 = n/2 and <j) = or n), Eqs. (Q])- 
d2J reduce to 

A = ^-, (5) 
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Hence, the post-SN orbital eccentricity is approximately 
proportional to the amount of mass lost during the SN ex- 
plosion if the spread in the values of the second term be- 
tween the square brackets in Eq. (O is sufficiently small. 
The question therefore arises whether typical values of 
Aq/[G(M() — M2)] for DNS progenitors are small enough 
to relax the requirement for low kicks to explain the ob- 
served Psp m -e relation. 



POPULATION SYNTHESIS 



To examine the role of polar kicks in establishing the 
P S pi n -e relation, we use the StarTrack binary population 
synthesis code to construct populations of DNSs for ming 
through the standard evolutionary channel (see ll24l l25ll 
for an extensive and detailed description of the popula- 
tion synthesis code). In this analysis, we consider popula- 



tion synthesis models with the same input parameters as 
the standard model (model A) of Belczynski et al. B25I1 . 
except that we do not allow for hypercritical accretion 
onto the NS during common envelope evolution and we 
modify the adopted kick velocity distribution to explore 
the effects of polar kicks. In particular, the first-born NS 
is assumed to be born with a kick velocity typical of 
isolated radio pulsars drawn from a Maxwellian veloc- 
ity distribution with dispersion a = 300kms _1 ll26[l27ll . 
while the second-born NS is assumed to be born with 
a kick velocity drawn from a Maxwellian velocity dis- 
tribution with dispersion a = 20km s (model K20) or 
a = 30010ns- 1 (model K300). The direction of the kick 
imparted to the first-born NS is furthermore assumed to 
be isotropic, while the direction of the kick imparted to 
the second-born NS is constrained to be within 10° from 
the progenitor's rotation axis. 

In the context of the standard evolutionary channel, 
the first-born NS is spun-up by dynamically stable mass 
transfer from the stripped down helium core of the NS's 
companion. The accretion rate is at all times limited to 
the Eddington rate for helium-rich matter. Spin-up is then 
modeled assuming the NS has an initial spin period of 1 s 
at the start of the mass transfer phase and accretion of a 
mass AM 1 increases NS's spin angular momentum by 



A/ S pin = ^GMi R AA AM U 
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wher e Ra j is the Alfven radius of the first-born NS (see 
also Il20l0 . We calculate the latter by randomly drawing 
a magnetic field strength from a uniform distribution 
between 10 9 and 10 10 G. 

A scatterplot of the present-day orbital eccentricities 
and spin-periods of the first-born NS in the resulting 
DNS populations is shown in Fig. [3] In the model where 
kicks typical of isolated radio pulsars (model K300) are 
imparted to the second-born NS, no correlation appears 
between P sp j n and e. However, when low kicks are im- 
parted to the second-born NS (model K20), the low-e 
end of the correlation is well sampled by the simulations, 
although a significant number of systems is also formed 
with spin periods P sp j n < 20 s below those of the observed 
systems. A possible cause may be the oversimplified na- 
ture of our NS spin-up model. These conclusions are in 
agreement with those of Dewi et al. 1 20] who studied the 
P sp i n -e correlation assuming isotropic kick distributions. 

As can be seen from Fig. [3] the main problem of 
the low kick scenario is its inability to explain the two 
observed DNS binaries PSRB 19 13+16 and PSRJ1811- 
1736 with eccentricities e > 0.6. The problem arises from 
the low NS progenitor masses (Mo < 2.7M ) at the end 
of the mass transfer phase spinning up the first NS: for 
Mq < 2.7 M and M\ = M2 = 1 .35M , the maximum at- 
tainable orbital eccentricity without a SN kick is 0.5. In 
the absence of kicks, high-e systems are therefore natu- 
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FIGURE 3. Scatterplots showing the present-day orbital eccentricities and spin-periods of the first-born NS in DNS binaries for 
the considered SN kick models (left: model K300; right: model K20). The present-day orbital eccentricities are determined by 
evolving the post-SN orbital elements under the influence gravitational wave emission up to the present epoch. Squares indicate the 
orbital eccentricities and spin-periods of the first-born NSs of observed Galactic DNSs. 



rally filtered out by the standard evolutionary channefl 
However, higher-mass progenitors for the second-born 
NS are possible if spin-up of the first-born NS can oc- 
cur during the common envelope phase of its compan- 
ion. This is illustrated in Fig.|4]where we show a scatter- 
plot of the present-day orbital eccentricities and progen- 
itor masses of the second-born NS in DNS binaries for 
the low kick model K20. In addition to the DNSs form- 
ing through the standard evolutionary channel, the fig- 
ure also shows those that form without any mass transfer 
following the common envelope phase. Since we do not 
consider hypercritical accretion, the first-born NS in the 
latter group of systems is not spun-up in our models. 

Figure |4] shows a clear dichotomy between low-e 
systems associated with low-mass progenitors for the 
second-born NS, and high-e systems associated with 
higher-mass progenitors for the second-born NS. The 
systems with low-mass progenitors all undergo stable 
mass transfer when the helium star evolves from core to 
shell helium burning and swells up to giant dimensions. 
At the start of the mass transfer phase these helium stars 
typically have masses below 3.5M . Stars more massive 
than 3 .5 Mq do not expand much after exhausting their 
central helium supply, thereby avoiding any further mass 



1 Dewi et al. [20] found a significant number of high-eccentricity 
systems in their population synthesis study of DNS binaries based on 
isotropic SN kicks with low kick velocities. While it is tempting to 
ascribe this difference to the different formation channels considered 
in our and their paper, both channels rely on spin-up of the first-born 
NS during a dynamically stable mass transfer phase prior to the SN 
explosion forming the second NS. The origin of the high-eccentricity 
systems in the absence of significant SN kicks in the Dewi et al. 12011 
study is therefore unclear, unless the final donor star masses at the 
end of mass transfer in the latter authors' simulations are substantially 
higher than those in our simulations. 
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FIGURE 4. Scatterplot showing the present-day orbital ec- 
centricities and progenitor masses of the second-born NS in 
DNS binaries for SN kick model K20. The arrows on the left 
indicate the orbital eccentricities of observed Galactic DNSs. 

transfer after the dynamically unstable common envelope 
phase. Since the eccentricities in the low kick model are 
predominantly determined by the mass lost during the 
SN explosion forming the second NS (and the subse- 
quent evolution towards lower eccentricities due to grav- 
itational wave emission), the separation of the second 
NS's progenitors in a low- and high-mass group causes in 
a dichotomy in the DNS orbital eccentricity range. While 
the observed sample of DNSs is still small, it is interest- 
ing to note that the observed DNS eccentricies show a 
lack of systems with 0.3 < e < 0.6. 

Our limited understanding of NS accretion during 
common envelope evolution prohibits us from extend- 
ing the simple spin-up model adopted in the case of dy- 
namically stable mass transfer to dynamically unstable 
mass transfer. While it has been suggested that NSs in 



a common envelope may accrete enough matter to col- 
lapse into a black hole J28ll29l[30l[3ll[32ll . counter argu- 
ments claim rotation and jets may prevent significant ac- 
cretion of matter 133113411 . Until significant improvements 
are made in our understanding of common envelope evo- 
lution, we are therefore left to speculate on the spin-up 
of the first-born NSs in the high-e systems in Fig. [4] In 
this context, we note that if NSs are able to accrete dur- 
ing common envelope phases and if the amount of matter 
accreted decreases with increasing mass of the donor star 
(cf. Fig |2), the group of high-e systems shown in Fig. [4] 
can be expected to form an extension of the P S pin _ e cor- 
relation found for low-e systems in the right-hand panel 
of Fig. [3] The presence of high-e systems in the observed 
P sp i n -e relation can be interpreted as support for a com- 
mon envelope model satisfying these assumptions. 

DISCUSSION AND CONCLUSIONS 

We investigated the role of spin-kick alignment in es- 
tablishing a correlation between the spin period of the 
first-born NS and the orbital eccentricity of DNSs us- 
ing the StarTrack binary population synthesis code and 
a simple prescription for the spin-up of a NS due to 
mass accretion. For DNSs forming through the standard 
evolutionary channel, spin-kick alignment is compatible 
with the observed P sp i n -e relation, but does not allevi- 
ate the requirement for low kick velocities proposed in 
previous population synthesis studies based on isotropic 
kick distributions. This is puzzling considering the strin- 
gent lower limits on the kick velocities imparted to the 
second-born NS in PSRB1534+12 and PSRB1913+16. 

Moreover, if low kicks are imparted to the second- 
born NS in DNSs, the standard formation channel can- 
not explain the formation of the high-eccentricity sys- 
tems PSRB1913+16 and PSRJ181 1-1736. A possible 
resolution would be a dichotomous formation channel 
where low-eccentricty DNSs are formed through the 
standard formation channel with low kicks imparted to 
the second-born NS, and high-eccentricity systems are 
formed through a formation channel where no mass 
transfer occurs after the common envelope phase of the 
second NS's progenitor and "normal" kicks typical of 
isolated radio pulsars are imparted to the second-born 
NS. The low-eccentricity systems on the P sp i n -e relation 
can be formed with either polar or isotropic kicks. How- 
ever, to obtain a P S p m -e relation at high eccentricities, 
some spin-kick alignment is required to counteract the 
increased spread in the post-SN orbital eccentricities in- 
troduced by the larger kicks. We will explore this in the 
continuation of this investigation. 
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